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Background:Experiments designed to study rare processes, such as neutrinoless double beta 
decay (0^/3/3), are crucial tests for physics beyond the standard model. These experiments rely on 
reducing the intrinsic radioactive background to unprecedented levels, while adequately shielding 
the detectors from external sources of radioactivity. Purpose: An understanding of the potential 
for neutron excitation of the shielding and detector materials is important for obtaining this level 
of sensitivity. Methods: Using the broad-spectrum neutron beam at LANSCE, we have measured 
inelastic neutron scattering on ™ at Cu. The goal of this work is focused on understanding the back- 
ground rates from neutrons interacting in these materials in regions around the Q-values of many 
candidate Ov/3/3 decay isotopes, as well as providing data for benchmarking Monte Carlo simula- 
tions of background events. Results: We extracted the level cross sections from the 7 production 
cross section for 46 energy levels in nat Cu . These level cross sections were compared with the 
available experimental data, as well as the ENDF/B-VII evaluation for discrete levels. Conclu- 
sions: For energy levels above 2 MeV we found significant discrepancies between the suggested 
level cross sections for both nuclei and our data. We found reasonable agreement between our 
measurement and the ENDF/B-VII evaluation for the total neutron inelastic cross section in 63 Cu. 
Our measurement of the total neutron inelastic scattering cross section in 65 Cu was 30% lower than 
the ENDF/B-VII evaluations, which we attribute to unobserved transitions in 65 Cu. Furthermore, 
we found that the implementation of the ENDF/B-VII evaluation in simulations did not properly 
model the decay properties of the nucleus to the degree necessary for estimating backgrounds in 
rear-event searches. Finally, we examined the potential implications of our measurements on Ov/3/3 
measurements and found that many of the commonly studied Ou/3/3 isotopes had Q-values below 
the cutoff for ENDF/B-VII evaluated discrete levels in either Cu nucleus. 



PACS numbers: 23.40.-S,25.40.Fq 

I. INTRODUCTION 

Neutrinoless double-beta decay (Qv(3/3) plays a key role 
in understanding the neutrino's absolute mass scale and 
particle-antiparticle nature [U? -[5] . If this nuclear decay 
process exists, one would observe a mono-energetic line 
originating from a material containing an isotope subject 
to this decay mode. One such isotope that may undergo 
this decay is 76 Ge. An experiment using germanium- 
diode detectors fabricated from material enriched in 76 Ge 
was the first experiment to established a half-life limit 
and a restrictive constraints on the effective Majorana 
mass for the neutrino [6l [7]. One analysis [8] of the 
data in Ref. [7] claims evidence for the decay with a 
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half-life of 2.23tg;^ x 10 25 y. Planned Ge-based 0vf3f3 
experiments [HI HD] will test this claim. Eventually, these 
future experiments target a sensitivity of >10 27 y or 
~1 event /ton- year to explore neutrino mass values near 
that indicated by the atmospheric neutrino oscillation 
results [IT]. 

The key to these experiments lies in the ability to re- 
duce intrinsic radioactive background to unprecedented 
levels and to adequately shield the detectors from ex- 
ternal sources of radioactivity. Previous experiments' 
limiting backgrounds have been trace levels of natural 
decay chain isotopes within the detector and shielding 
components. The 7-ray emissions from these isotopes 
can deposit energy in the Ge detectors producing a con- 
tinuum, which may overwhelm the potential Qvj3j3 sig- 
nal peak at 2039 keV. Great progress has been made 
identifying the location and origin of this contamination, 
and future efforts will substantially reduce this contribu- 
tion to the background. The background level goal of 
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1 event/ton-year, however, is an ambitious factor better 
than the currently best achieved background level for Ge- 
based experiments [7J . If the efforts to reduce the natural 
decay chain isotopes are successful, previously unimpor- 
tant components of the background must be understood 
and eliminated. The work of Mei and Himejl2] recog- 
nized that (n, n 1 "/) reactions will become important for 
ton-scale double-beta decay experiments. 

Reference [13] recognized that the specific 7 rays 
from Pb isotopes at 2041 and 3062 keV are particu- 
larly troublesome. The former is dangerously near the 
2039.00±0.05-keV Q-value for zero-neutrino double-beta 
decay in 76 Ge and the latter can produce a double-escape 
peak line at 2040 keV. That paper pointed out that the 
cross sections to produce these lines in nat Pb were un- 
measured and hence set to zero in the data bases of the 
simulation codes used to design and analyze 0i/(3/3 data. 
This result indicated the importance of assessing all ma- 
terials used in Op/3/3 experiments for such problematic 
lines. In this work we examine Cu in this context, be- 
cause it is a material frequently used in large quantities 
in Ov/3/3 experiments. 

Previous authors have studied (n,n'7) reactions in Cu 
usually using both natural Q3HIE] [EH2S] and isotopic 
samples [27 32 to isolate the isotopic effects. Previous 
measurements, however, have not extended to nuclear ex- 
citation energies greater than 1.5 MeV and therefore the 
production rates of 7 rays beyond this energy were not 
measured. This paper presents measurements of nat Cu 
(n,n'7) production cross sections of the 7 rays for nat Cu. 
Although our work was motivated by neutron reaction 
considerations in materials that play important roles in 
the Majorana [33J design, the results have wider utility 
since Cu is used by numerous low-background experi- 
ments. 



II. EXPERIMENTAL PROCEDURE AND DATA 
ANALYSIS 

The experimental data were collected at the Los 
Alamos Neutron Science Center (LANSCE) located at 
Los Alamos National Laboratory. Spallation neutrons 
are produced from an 800-MeV proton beam interact- 
ing with a nat W spallation target in the Weapons Neu- 
tron Research (WNR) Facility [33]. The incident proton 
beam is delivered in short pulses spaced 1.8 /zs for 625 
/is at a rate of 40 Hz. The resulting neutrons range in 
energy from 0.2-800 MeV. To ensure that the beam area 
illuminates only the Cu target, the beam is collimated 
to a 1.9-cm radius with a lead collimator. The neutron 
flux is monitored using an in-beam fission chamber with 

235,238tj foUs 

The ™ at Cu target was placed at the center of GEANIE 
spectrometer [35], 20.34-m downstream from the produc- 
tion spallation target, and 60° to the right of the proton 
beam. Three stacked foils of na< Cu were placed perpen- 
dicular to the neutron beam. The foils, each measur- 



ing 50 x 50 x 0.5 mm, had a total mass of 139 g. The 
surrounding GEANIE spectrometer comprises 26 germa- 
nium detectors; 16 coaxial detectors with an energy range 
of four MeV, and 10 Compton suppressed planar detec- 
tors with an energy range of one MeV. During the run 
cycle, the detectors experience significant neutron expo- 
sure, and subsequently some of the detectors have re- 
duced energy resolution, or gain instabilities. During the 
course of the experiment three planar detectors and six 
coaxial detectors (nine total detectors) were not used for 
exhibiting such problematic behavior. The remaining six 
planar and ten coaxial detectors were used for this anal- 
ysis. 

The procedure for extracting 7-ray production cross 
sections from the GEANIE spectrometer has been dis- 
cussed in great detail in several papers (see for example 
[33(37].). The GEANIE spectrometer provides the yields 
and energies of discrete cascade 7 rays. Yields are con- 
verted into partial 7-production cross-sections using cali- 
brated detector efficiency, deadtime corrections, neutron 
flux, and target thickness. 



III. RESULTS AND DISCUSSION 

The 7 production cross sections were the primary mea- 
sured quantities of the present experiment. In total 119 
7-ray transitions were measured in 63 > 64 > 65 Cu for neu- 
tron energies from 0.2-100 MeV. The total inelastic and 
the level cross sections for individual excited states were 
deduced from the evaluated level scheme of 63 Cu, 64 Cu, 
and 65 Cu. The level cross sections for selected levels in 
both Cu nuclei are presented in the sections below, as well 
as the total neutron inelastic cross sections and (n,2n7) 
cross sections in Cu. Since the goal of this work is to 
place these measurements in the context of the typical 
background simulation for a neutrinoless double beta de- 
cay experiment, we begin our analysis by presenting our 
data in the context of data from the literature, and the 
evaluated cross sections presented in the ENDF/B-VII 
evaluation on Cu. 

The ENDF/B-VII presents cross sections for the exci- 
tation of an energy level for a range of incident neutrons. 
At high incident neutron energies discrete gamma rays 
may be produced by various excitation paths, particu- 
larly feeding from 7 decay of higher lying levels. The 
neutron inelastic partial cross sections to discrete lev- 
els supplied by ENDF/B-VII do not include this impor- 
tant component. These level excitation cross sections 
differ from those obtained through individual 7 produc- 
tion cross sections for a given excited state because the 
7 production cross sections include the effects of feed- 
ing from higher energy levels. Thus level cross sections 
deduced from summing up the individual 7 production 
cross sections associated with a given level are a measure 
of the level de-excitation cross section. We have chosen 
to present these two types of level cross sections together 
because the focus of this paper is the accurate simulation 
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of both the excitation and the de-excitation of the energy 
levels in Cu. 

The second section of the results discusses the results 
of this experiment in the context of two common statis- 
tical model calculations, TALYS[3E] & CoH3[3H]. Statis- 
tical model calculations are quite important to the data 
evaluation process, as they serve to fill in the gaps of ex- 
perimental data and provide cross sections where no mea- 
surements are available. Thereafter we compare the nat- 
ural Cu(n,n'7) and (n,2n7) cross sections with GEANT4 
simulations. Finally we discuss the implication of this 
measurement on neutrinoless double beta decay measure- 
ments. 



Comparison with data and the ENDF/B-VII 
evaluations 



We extracted the level de-excitation cross sections from 
the 7 production cross section for 46 energy levels in 
nat Cu . These level de-excitation cross sections were com- 
pared with the available experimental data, as well as the 
ENDF/B-VII evaluation for level excitation cross section 
of discrete levels. Figures [2j and [3] show the level de- 
excitation cross sections for the first energy levels in 63 Cu 
and 65 Cu. In these Figs, we show our data along side the 
ENDF/B-VII evaluations for these levels, as well as sev- 
eral (n,n'7) measurements (circles) [Tj2 [22 [23 US] . The 
second energy levels in both nuclei have very weak decay 
branches to the first excited level, therefore the strength 
from these indirect excitations doesn't begin to influence 
the level de-excitation cross section until the neutron en- 
ergy exceeds the energy required to excite the nucleus 
into the third excited level (E x =1326 keV for 63 Cu and 
1481 keV for 65 Cu, see Fig. [l] for the level diagram of 
both Cu nuclei). The level de-excitation cross sections 
peak near 1.34, and 1.689 MeV for 63 Cu and 65 Cu, re- 
spectively. For neutron energies below the maximum 
cross section there is minimal contribution from these 
higher energy levels, and therefore the level de-excitation 
cross sections can be entirely attributed to direct excita- 
tion. The measured maximum level de-excitation cross 
sections are consistent with the ENDF/B-VII evaluation. 
The maximum level de-excitation cross section for 63 Cu 
is o" ma2; =226(14) mb, which is roughly 6% higher than 
a max =212{8) mb for 65 Cu. 

Fig. [4] shows our experimental data for the fourth and 
fifth energy levels in both nuclei, together with previ- 
ous measurements, as well as the ENDF/B-VII evalu- 
ation for the four levels shown. The cross sections for 
the energy levels are determined by summing the 7 pro- 
duction cross sections for each of the decay branches. 
In the case of the 1326-keV level of 63 Cu, the level de- 
excitation cross section is determined by summing the 7 
production cross sections for the 1412, 742, and 442-keV 
7-ray transitions. The level de-excitation cross section 
for the 1623-keV level in 65 Cu was obtained by summing 
the 7 production cross sections of the 1623, 852.7, and 
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FIG. 1. The level scheme of the six lowest energy levels in 
63 Cu ED] and 65 Cu ST]. 
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FIG. 2. The total 7-ray production cross section for the first 
excited level (E;c=669 keV) in 63 Cu from this measurement, 
and several (n,n'7) measurements (squares) [I9l l22l [23l I26] . 
the ENDF/B-VII evaluation (solid line), the TALYS1.2 and 
CoH3 calculations (dashed lines) for this level. At high in- 
cident neutron energies discrete 7 rays may be produced by 
various excitation paths, especially feeding from 7 decay of 
higher lying levels, and the neutron inelastic partial cross sec- 
tions to discrete levels. Note that the ENDF/B-VII level cross 
sections do not include feeding from higher lying levels that 
is in the measured cross sections and nuclear model calcu- 
lations. The level de-excitation cross sections presented for 
the statistical models are similar to those presented from this 
measurement; they are calculated by summing up the 7 pro- 
duction cross section for a particular level. 
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FIG. 3. The total 7-ray production cross section for the first 
excited level (E;c=770 keV) in 65 Cu from this measurement, 
as well as several (n,n'7) measurements (squares) |22l I26] , 
the ENDF/B-VII evaluation (solid line), the TALYS1.2 and 
CoH3 calculations (dashed lines) for this level [42]. For a 
description of the types of level de-excitation cross sections 
presented here please see Figure [2] 



507.9 keV transitions. The 1547-keV level de-excitation 
cross section was obtained by summing the 7 production 
cross sections for the 1547, 877, and 584-keV 7-ray tran- 
sitions. Finally the 1725-keV level de-excitation cross 
section in 65 Cu was determined by summing the 7 pro- 
duction cross sections for 1725.92, 954.5, and 609.5-keV 
7-ray transitions. The branching ratios for each of these 
energy levels was calculated from the ratio of the 7 pro- 
duction cross section for each transition to the total level 
de-excitation cross section. In all cases the branching ra- 
tios for the transitions were compared with those given 
in ENDSF and were generally found to be in agreement 
with the literature. Branching ratios were computed for 
all 81 neutron energy bins used in this analysis, and were 
constant across the entire neutron energy range. In gen- 
eral, the measured level de-excitation cross sections were 
in good agreement with previous experimental data, and 
for a limited energy-range with the ENDF/B-VII evalu- 
ation. 

Guenther et. al measured the 7 production cross sec- 
tions at two angles (55° and 125°) with incident neutrons 
up to 4.5 MeV, and [55]. We find that our measure- 
ments, which include a much larger solid angle, are about 
10% higher, although both measurements have compa- 
rable shapes. Similar measurements for the first excited 
level in 65 Cu do not show such a disagreement. Our mea- 



surements are consistent with 7 production cross section 
measurements from Nishimura et. al 22\ at the JAERI 
5.5 MV Van De Graaff, as well as the measurements of 
Tucker et. al [23 for 63 Cu. 

Figure [4] shows the measured level de-excitation cross 
sections along with the available data, and the ENDF /B- 
VII evaluated level excitation cross sections for the fourth 
and fifth energy levels in both nuclei. Our data agree with 
the limited data from measurements. The lack of exper- 
iment data is important because these partial neutron 
cross sections typically provide the data for the ENDF 
evaluations that give the discrete level excitation cross 
sections. Due to the difficulty of resolving closely spaced 
levels, neutron measurements, are usually available at 
only low energy levels and for a limited range of incident 
neutron energies. Prior to our work, the ENDF/B-VII 
evaluation on Cu used very little experimental data above 
the fifth energy levels in either nuclei, and only included 
data for three incident neutron energies. 

Figs [5J [7] are a comparison of our measured level de- 
excitation cross sections in this experiment with the level 
excitation cross sections recommended by ENDF/B-VII 
where there was no experimental data available at the 
time of the evaluation. Feeding from 7-rays at higher 
levels tends to have pronounced effects around 4-5 MeV, 
well above the peak cross section for the states discussed 
in this paper. In the energy region at or below the maxi- 
mum cross section, the population mechanism for the lev- 
els is assumed to be direct excitation and therefore is di- 
rectly comparable with the ENDF/B-VII recommended 
level excitation cross sections. Comparing the maximum 
cross sections of these levels with our experimental data 
reveals the issue with using the evaluation at higher en- 
ergies; the ENDF/B-VII evaluated discrete level cross 
sections at higher energy are typically determined by sta- 
tistical models, and at least in this evaluation many of 
the evaluated discrete level cross sections disagreed with 
the experimental cross sections by more than 10%. 

Of the three levels presented here only one of the lev- 
els was in agreement with the ENDF/B-VII evaluation. 
The decay of these energy levels produces 7 rays that 
are in the vicinity of the expected neutrinoless double 
beta decay signal for several neutrinoless double beta 
decay isotopes. In the case of 2081-keV level in 63 Cu, 
the ENDF/B-VII evaluation overestimates the maximum 
cross section by about 25%. This level shows significant 
strength at the higher energies from feeding from higher- 
lying levels. The evaluation, however, does not directly 
include the feeding effects for this level because the evalu- 
ation does not include discrete higher lying levels that di- 
rectly feed this level. In most simulation packages invok- 
ing the ENDF/B-VII data set, the feeding of this level is 
handled predominantly by the statistical models. These 
models typically don't cascade through the defined levels, 
but rather through a statistically determined level den- 
sity. The 2535-keV level shows much better agreement 
with the ENDF/B-VII evaluation, the maximum cross 
sections are consistent with the measured data, and be- 
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FIG. 4. The total 7-ray production cross section for the fourth (E K =1412 & 1623 keV) and fifth (E x =1547 & 1725 keV) energy 
levels in 63 Cu (left column) and 65 Cu (right column) from this work as well as the inelastic neutron scattering cross sections 
from [T51[21], and the ENDF/B-VII evaluation for these levels [32]. The level de-excitation cross sections from the TALYS1.2 
and CoH3 calculations are presented as well (dashed lines) for these level. For a description of the types of level de-excitation 
cross sections presented here please see Figure [2] 



cause there is not significant feeding from higher lying 
states the shape of the level cross sections are in agree- 
ment with each other. The maximum cross section of the 
2808-keV level in 63 Cu, however is about 50% lower than 
the ENDF/B-VII evaluation. Overall, our experimental 
data found disagreement in 13 of the 24 ENDF/B-VII 
evaluated levels above 2 MeV: eight had maximum cross 
sections lower than the ENDF/B-VII estimate, three had 
maximum cross sections values higher than the evalua- 
tion, and two of the levels were not seen. In all cases the 



discrepancies were larger than 10%, and in one particular 
case, the ENDF/B-VII the maximum cross section was 
overestimated by a factor of three. 

In addition to comparing our measured level de- 
excitation cross sections with the ENDF/B-VII evaluated 
level excitation cross sections, we also compared our data 
with the ENDF/B-VII evaluated total neutron inelastic 
cross sections for Cu and 65 Cu. The total neutron in- 
elastic cross sections were obtained by summing the 7 
production cross sections for the ground state transitions 
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FIG. 5. The total 7-ray production cross section for the 
2081-keV level in 63 Cu from this measurement, as well as the 
ENDF/B-VII evaluation for this excited level. For a descrip- 
tion of the types of level de-excitation cross sections presented 
here please see Figure [2] The decay of this level produces a 7 
ray in the vicinity of the 76 Ge endpoint. 



FIG. 6. The total 7-ray production cross section for the 
2535-keV level in 63 Cu from this measurement, as well as the 
ENDF/B-VII evaluation for this level. For a description of 
the types of level de-excitation cross sections presented here 
please see Figure [2] The 2535-keV level is in the vicinity of 
the 130 Te endpoint energy. 



in each isotope, and are shown in Figs. [H]and[9j The total 
inelastic cross sections deduced from 7 production cross 
sections are comparable to the ENDF/B-VII evaluation 
provided there isn't significant strength associated with 
the ground state decays of energy levels greater than 4 
MeV. In the case of 63 Cu, the total neutron inelastic cross 
section is in excellent agreement with both experimental 
data [201 HQ H3 133], and the ENDF/B-VII evaluation. 
Around 20 MeV there is a slight increase in total neu- 
tron inelastic cross section from this measurement which 
is due to the 65 Cu (n,3n7) reaction. According to the 
description of the ENDF/B-VII evaluation, 26 levels be- 
low 3.7 MeV were used to build the level information for 
63 Cu 44J. The levels used in the total neutron inelas- 
tic cross sections are the same below 2.889 MeV. The 
ENDF/B-VII evaluation, however includes three levels 
at 3.3, 3.48, and 3.7 which are assumed to be collective 
excitations and provide significant contributions to the 
inelastic-scattering and 7-ray production cross sections 
[44] . There was no evidence for any ground-state transi- 
tions from these levels, nor transitions to any of the levels 
assumed in the decay paths provided in the ENDF/B-VII 
evaluation. It is difficult to evaluate the existence of these 
levels, since the exact excitation energies, spins, parities 
and decay paths are unknown. There have been sev- 
eral inelastic proton experiments that have seen levels at 
these energies with quite considerable strength [351 ES] > 



however there are no 7-ray transitions in the current ex- 
perimental data with sufficient strength to indicate the 
presence of these levels. 

The total neutron inelastic cross section for 65 Cu (see 
Fig. [9| is considerably lower than the ENDF/B-VII eval- 
uation. For this evaluation, the 7 production cross sec- 
tion for 19 levels below 3.360 MeV were included in the 
summation. As with the 63 Cu evaluation, four collec- 
tive levels were included in the ENDF/B-VII evaluation 
at 2533, 3080, 3350, and 3500 keV. In addition to these 
three collective levels, a forth presumed collective level 
at 2533 keV was also not observed in our work. The 
ground-state transitions were also not observed for these 
levels, nor was there any enhanced 7-ray transition that 
would indicate the existence of these strong collective ex- 
citations. The total neutron inelastic cross section from 
this experiment is in agreement with the measurements 
of [20 , but are somewhat lower than the results of [25], 
and [23] . In general, the total neutron inelastic cross sec- 
tion for 65 Cu is about 30% lower than the ENDF/B-VII 
evaluation, and the [231 [35] measurements. 

From the current measurement, there appears to be a 
significant discrepancy between the high-energy states in 
these two nuclei. There are several potential explanations 
for the deficit in the continuum strength of 65 Cu It is 
quite likely that the experimental value is low because 
of unobserved gamma transitions. Either the gammas 
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FIG. 7. The total 7-ray production cross section for the 
2808-keV level in 63 Cu from this measurement, as well as the 
ENDF/B-VII evaluation for this level. For a description of 
the types of level de-excitation cross sections presented here 
please see Figure [2] The 2808-keV level is in the vicinity of 
the 116 Cd endpoint energy. The range of discrete levels in the 
Talysl.2 simulation did not extend high enough in energy to 
provide a level de-excitation cross section for this state. 
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FIG. 8. The total neutron inelastic cross section for Cu, as 
well as the results of (n,n') measurements [43] (squares), the 
(n,n'7) measurement (circles), and the ENDF/B-VII 

evaluation for the total neutron inelastic cross section. 



FIG. 9. The total neutron inelastic cross section for Cu, as 
well as the results of an (n,n') measurements [2"4] (squares), the 
(n,n'7) measurements [13 HI] (circles), and the ENDF/B-VII 
evaluation for the total neutron inelastic cross section. 



are at higher energies, above the 4 MeV cutoff of the 
GEANIE spectra, or the decay is fragmented into many 
weaker gammas below the detection limits, or both. A 
long-lived isomer in 65 Cu could also remove some of the 
cross section, because GEANIE measures mostly prompt 
gammas, but no such isomers are known in 65 Cu. The 
most likely explanation for the decrease in the continuum 
strength is unobserved 7 transitions. Such unobserved 
strength is common in odd-odd nuclei close to the mid- 
shell; nuclei near the mid-shell region typically experience 
maximal deformation |47j . 



B. Comparison with Statistical Models 

Two different codes were used for the statistical 
model calculations in this analysis, TALYS-1.2|38j. and 
COH3 39J. These codes are both based on the Hauser- 
Feshback formalism, and both rely on local optical model 
parameters for the calculations. The statistical models 
were used to produce 7 production cross sections for sev- 
eral transitions in 63 Cu and 65 Cu. These 7 production 
cross sections were summed to give level de-excitation 
cross sections similar to those presented throughout the 
text. The results of these calculations were compared 
with the first, fourth and fifth energy levels in both 63 Cu 
and 65 Cu, as well as the 2081-keV level in 63 Cu. 

Examining the agreement between the various mod- 
els and the experimental data, it becomes very apparent 



that neither model is able to predict the strength of in- 
dividual levels to within 10% of the experimental data. 
For the first energy levels in both Cu nuclei, both models 
provide a reasonable approximation of the actual cross 
section up to the maximum cross section. The CoH3 cal- 
culations overpredict the effects of feeding from higher 
lying energy levels to the first excited level. In fact for 
63 Cu the CoH3 overpredicts the cross section by about 
20% between 5 and 15 MeV, while the TALYS-1.2 model 
underpredicts the cross section by about 30% in the same 
energy region. While for 65 Cu, both models overpredict 
the cross section in this region by about 55%. Overall, 
both codes overpredict the summed strength in 65 Cu by 
about 50%, while for 63 Cu the summed strength for the 
CoH3 model is about 11% higher than our data, and the 
TALYS-1.2 calculation is about 25% lower than our data. 

At higher energy levels there are similar problems. 
Fig. [4] shows the third and fourth excited levels for 
both Cu nuclei. The maximum cross sections for the 
third and fourth excited levels are overpredicted by both 
calculations for both nuclei. The fourth and fifth ex- 
cited levels in 63 Cu are overpredicted by 6-8%, and 14- 
17%, respectively. While the same levels in 65 Cu are 
overpredicted by 25% in both calculations. The CoH3 
model achieves reasonable agreement with the fourth 
excited level in 63 Cu, but has 67% less total strength 
than our experimental data for the third excited level. 
The TALYS-1.2 calculation underpredicts the summed 
strength in both excited levels: having roughly 90% of 
the observed total strength in the third excited level, 
and 85% of the observed total strength in the fourth ex- 
cited level. The third excited level in 65 Cu is in reason- 
able agreement with both models; the CoH3 calculation 
slightly underpredicts the summed strength by about 
10%, while the TALYS-1.2 calculation overpredicts the 
summed strength by 11%. Neither calculation is in agree- 
ment with the fourth excited level in 65 Cu; both overes- 
timate the summed strength by 260-290%. Finally, we 
turn to the 2081-keV level in 63 Cu. This level is impor- 
tant for experiments like Majorana because the pri- 
mary transition emits a 7-ray in the vicinity of the 76 Ge 
endpoint energy. For this level the calculations overpre- 
dict the maximum cross section by 6-8%. The TALYS-1.2 
calculation underpredicts the strength in the continuum 
region (taken to be between 5-15 MeV) by about 6%, 
while the CoH3 calculation is roughly in agreement with 
the experimental data. 

The statistical models typically in use for estimating 
7 production cross sections for individual excited levels 
display inconsistencies from nucleus to nucleus, and even 
among excited levels in the same nucleus. Neither model 
seems to be able to accurately predict the strength in a 
particular level to the precision necessary for double-beta 
decay and other high-sensitivity, low-background exper- 
iments. The models vary considerably in their model- 
ing of a particular level; showing agreement to within 
10% in one level and then over predicting the summed 
strength in a neighboring level by a factor of 3. In ad- 
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FIG. 10. The total inelastic cross section for nat Cu (circles), 
the ENDF/B-VII evaluation (line), and the results of the 
GEANT4 simulation (squares). 



dition, both codes seem to overpredict the contributions 
from the continuum in 65 Cu. Experiments designed to 
search for rare signals require a better understanding of 
their background contributions than can be afforded with 
the current statistical model calculations. 



C. Monte Carlo Simulations 

For most f3/3 studies the successful observation of a 
0u/3/3 signal relies on the ability to achieve an extremely 
low-background experimental setup. As this is the first 
generation of experiments to attempt such a low sensi- 
tivity, the estimates of the background levels are highly 
dependent on simulations. To quantify the accuracy 
of these simulations, we compared our experimentally 
determined total inelastic cross sections with those ex- 
tracted from simulations. For our analysis we simulated 
ten-million neutrons scattering off a 139-g block of no *Cu 
with the same dimensions as our actual target. The neu- 
trons had a flat energy distribution between to 40 MeV. 

The simulations were done in GEANT4-09-03 using 
the QGSPJLEJBICTSO physics list for inelastic neutron 
scattering [48j? ]. The ENDF/B-VII libraries were called 
to account for both elastic and inelastic neutron scatter- 
ing below 20 MeV. Above 20 MeV, the elastic scattering 
process was modeled using the low-energy elastic model, 
while the inelastic process was modeled from 19.9 MeV to 
9.9 GeV with the Binary Cascade model [IS]. The results 



of the simulation are shown in Figs. 10 and 11 The curve 



o 

<D 

CO 



O 

O 

as 
o 



a, Cu(n, 2 n y) E x < 644 keV 

• Present Data 
■ QGSP_LE_BIC_ISO 
ENDF/B-VII 




10 20 30 

Neutron Energy [MeV] 

FIG. 11. The (n,2n) cross section for nat Cu(circles), 
the ENDF/B-VII evaluation (line), and the results of the 
GEANT4 simulation (squares). 
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FIG. 12. Comparison of spectra produced from GEANT4 
simulation with that obtained from the present experiment. 
The top figures shows the experimental spectrum from this 
experiment, while the bottom spectrum shows the results of 
a GEANT4 simulation. The specifics of the GEANT4 simu- 
lation are outlined in the texts. 



marked as QGSP_LE_BIC_ISO refers to the summation 
of all inelastic scattering processes (n, n'7) that occur for 
a particular incident neutron, and normalized to the total 
number of simulated neutrons at that particular incident 
energy. The red circles, denoted as the present data, are 
the total inelastic cross sections for ™ at Cu as a function 
of incident neutron energy. These data points represent 
the summed strength of 25 levels and 14 levels below 3.1 
MeV in 63 Cu and 65 Cu, respectively. The ENDF/B-VII 
evaluation for nat Cu is depicted by the solid purple line, 
and is the sum of the natural abundance corrected 63 Cu 
and 65 Cu total inelastic cross sections. As can be seen 
from Figs. [lTJ| and [TT| GEANT4 simulations are consis- 
tent with the ENDF/B-VII evaluation. 

Geant4 implements the ENDF/B-VII evaluations in 
much the same way as other simulation programs. Once 
the inelastic process is activated a random number gener- 
ator is used to decide among the various available energy 
levels, which are weighted according to their contribu- 
tion to the total cross section. In the region where the 
ENDF/B-VII provides level excitation cross sections, and 
branching ratios, the simulations produce 7 rays that 
correspond with known levels and branching ratios in 
the particular nucleus. The issue comes when the ex- 
citation energy exceeds the ENDF/B-VII cutoff energy 
for discrete level, and the ENDF/B-VII continuum file 
is employed. The summed strength in this file usually 
accounts for upwards of 70% of the total strength in 
the ENDF/B-VII evaluation, and unlike the evaluated 
level data, doesn't have a decay path, rather it relies on 



model-generated 7 production probabilities, gamma-ray 
emission spectra, and neutron emission-energy spectra. 
At this point the decay no longer proceeds through de- 
fined levels in the nucleus, and the emitted 7 rays no 
longer correspond to transitions between nuclear levels. 
In addition given the significant strength associated with 
this file, this region becomes significant at very low inci- 
dent neutron energies; neutron above 3 MeV interacting 
in Cu are predominately described by the ENDF/B-VII 
continuum file. Figure [12] shows a comparison of the 
GEANT4 simulated spectrum compared with the exper- 
imental spectrum observed in the present experiment. 
The simulation shows significant strength in the 2 MeV 
region, which is not observed in the actual experimen- 
tal spectrum. The ENDF/B-VII evaluation appears to 
reproduce some 7-ray transitions very nicely. Although 
there appears to be some strange structure that appears 
below many peaks that is not seen in the actual experi- 
ment data (see for example the 2 and 2.1 MeV regions in 
Fig. ~ 



12) 



Furthermore, the simulation predicts strength 
from more states than are actually seen in the data, while 
also the relative strengths of several transitions do not 
appear to agree with the actual data. 



D. Cross Sections relevant to Ov/3/3 decay searches 

Cu is a very popular material in Qvj3fi decay searches. 
It has good thermal properties which allows its use in 
cooling systems. Cu can also be manufactured extremely 
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cleanly, and thus is generally used in large quantities 
around these low-background setups. One issue that has 
not been examined is the probability that neutrons inter- 
acting in Cu will produce 7 rays that might interfere with, 
or even replicate a Ov(3(3 signal. Only four of the nine 
potential Qv(3(3 isotopes listed in Table [I] had Q-values 
below the cutoff for ENDF/B-VII evaluated excited lev- 
els in either Cu nucleus. The strength associated with 
excited levels above the ENDF/B-VII cutoff energy have 
been folded into the continuum file, and thus a proper 
parameterization of the population, and decay of these 
excited levels will not be explicitly provided to simula- 
tions employing ENDF/B-VII. 

A quick literature search revealed several 7 ray lines 
that could potentially interfere with several isotopes fre- 
quently used in (3(3 studies (see Table [l}. In all, 12 7-ray 
lines were identified as potential background lines for the 
Ov(3(3 experiments. Of these 12 7-ray lines, seven were 
outside of the chosen energy range of the GEANIE de- 
tectors for this experiment. Three of these high-energy 
7-ray lines had detectable lower-energy 7-ray transitions 
which allowed for a cross section estimate. The energy re- 
gions surrounding the Q-value, Q-value + 511 keV (SEP, 
single-escape peak), and the Q-value + 1022 keV (DEP, 
double-escape peak) for the isotopes listed in Tab. [IJwere 
examined to rule out any additional unknown 7-ray lines 
that might not be currently noted in the literature. Ob- 
viously, the 7 ray lines with sufficient energy to produce 
DEP's near the Qpp of 82 Se, 100 Mo, and 150 Nd were too 
high in energy to be examined in this study. 



TABLE I. A list of frequently studied (3(3 isotopes and their 
Q-values. 7-production cross sections or upper limits for im- 
portant transitions in na *Cu are given for neutrons between 
2.897 and 4.196 MeV. Where the cross section is listed as NA, 
the experimenters were unable to place a limit due the 7-ray 
being outside the range of the detection system. 



(3(3 isotopes 


7 ray 


SEP 


DEP 




[mb] 


[mb] 


[mb] 


48 Ca 


65 Cu NA 


65 Cu NA 


65 Cu NA 


76 Ge 


65 Cu <0.388(3) NA 






82 Se 


63 Cu 9.42(32) 


<0.324(3) 


63 Cu NA 


96 Zr 


65 Cu 1.12(2) 


<0.241(3) 


65 Cu NA 


100 Mo 


63 Cu 9.42(32) 


63 Cu 0.59(22) 


63 Cu NA 


116 Cd 


63 Cu 4.41(23) 






iso Te 


B3 Cu 9.42(32) 


63 Cu 9.42(32) 


63 Cu <0.316(3) 


136 Xe 


63 Cu 0.62(10) 


< 0.392(3) 


63 Cu 1.03(10) 


150 Nd 


< 0.319(3) 


< 0.265(3) 


65 Cu NA 



a. 48 Ca Endpoint Energy The CANDLES experi- 
ment is an experiments looking to study Qv(3(3 in 48 Ca 
[50] , The endpoint energy for this isotope is above the 
energy range of the GEANIE spectrometer in this exper- 
iment, and thus cross section information is not available 
for this isotope. Recent measurements have shown that 
CaMo04 scintillators might be a viable option for study- 



ing 48 Ca; these types of detectors would potentially have 
a 4% energy resolution at 3-MeV when operated at 244 
°K [ST]. With this energy resolution, there are 35 7-ray 
transitions in Cu in a 340-keV energy region around the 
48 Ca endpoint energy. All of these 7-ray transitions de- 
rive from excited levels that are above the cutoff energy 
for the ENDF/B-VII evaluated discrete levels in either 
Cu nucleus. 

b. 76 Ge Endpoint Energy Experiments designed to 
study neutrinoless double beta decay in 76 Ge typically 
employ enriched Germanium diodes as an active detec- 
tion system. These crystals have excellent energy resolu- 
tion in the endpoint region - typically corresponding to a 
resolution of « 0.2% or about 4 keV [52]. The literature 
states that 65 Cu has a high energy level {E x = 3157(3) 
keV) that emits a 2041(3)-keV 7-ray line during the de- 
cay. The 2041-keV 7 ray line associated with the decay of 
this excited level was not observed in the GEANIE data, 
and consequently only a limit can be placed on the 7 pro- 
duction cross section. A review of the literature revealed 
that this is a weakly excited level, currently only detected 
in a limited number of experiments |29l I53j . Fortunately 
this decay is not to the ground state, and should be ac- 
companied by a 1115-keV 7-ray line. An examination 
of the 7 rays emitted in coincidence with the 1115-keV 
7-ray line did not yield any strength in the 2041-keV re- 
gion. Table [I] gives the sensitivity of our measurements 
to a 7 ray produced in this energy region. If this excited 
level exists and it decays via a 2041-keV 7 ray, it would 
likely have a maximum 7 production cross section below 
0.35 mb. There are no additional 7-ray transitions in 
the SEP & DEP energy region. 

c. 82 Se, 96 Zr, and 100 Mo Endpoint Energy The Se- 
lenium, Zirconium and Molybdenum isotopes have all 
been grouped together because experiments such as 
MOON and SuperNEMO are designed to measure any 
of these isotopes, in addition to several other isotopes in 
their setups [53] HSJ ■ There doesn't appear to be many 
interfering 7 rays lines in the immediate vicinity of 82 Se. 
There are several excited levels that emit 7 rays on either 
side of the 2995-keV endpoint energy in 82 Se, and the 7 
production cross sections are shown in Fig. [13] There 
is a single excited level that emits an interfering 7 ray 
within 5 keV of the 96 Zr endpoint energy, although this 
decay path is not the dominant decay channel for this 
energy level. The 7 production cross sections for this 
level, located at 3335 keV, is shown in Fig. 14 Both 
nuclei exhibit 7-ray lines that could double-escape into 
the respective regions of interest. Again in the case of Zr, 
the decay channel that emits a 4371-keV 7 ray from the 
7472.7-keV excited level is quite weak, occurring with a 
12% probability. 

In the immediate vicinity of the 100 Mo endpoint energy 
there is a ground-state transition from the 3032.70-keV 
energy level from 63 Cu. The 7 production cross section 
for this transition as a function of energy is shown in Fig. 
13 Expanding the search region to 3% energy resolution 
reveals an additional 14 7-ray transitions in raat Cu in the 
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FIG. 13. The 7 production cross sections for 2977, 3032, and 
3042-keV 7-rays in 63 Cu. These 7-ray are all ground state 
transitions, and therefore the energy of the 7-ray corresponds 
to the energy level in 63 Cu. These 7-ray transitions are in 
the vicinity of the 82 Se and 100 Mo endpoint energies. The 
corresponding levels for these 7-rays are not included in the 
ENDF/B-VII evaluation for 63 Cu. 



FIG. 14. The 7 production cross section for 3355-keV 7-ray 
in 65 Cu. This 7 ray transition is in the vicinity of the 96 Zr 
endpoint energy, and is a ground state transistion from the 
3355 energy level in 65 Cu. The corresponding level for this 
7-ray is not included in the ENDF/B-VII evaluation for 65 Cu. 



Mo region of interest. A similar search in the SEP and 
DEP regions revealed an additional 34 potential back- 
ground 7-ray lines. Again, as in the case of Ca, the 
excited levels that produce these 7-ray lines are above 
the cutoff energy for the ENDF-BVII evaluated discrete 
levels in either Cu nucleus. 

d. 116 Cd Endpoint Energy In the endpoint energy 
region of 116 Cd there is a ground state transition from an 
excited energy level in 63 Cu that is in the vicinity of the 
116 Cd endpoint energy. While the dominate decay path 
for this level is to the ground state, there are additional 
decay branches that might facilitate the identification of 
this as a potential background line. Fortunately, this 
level is included in the ENDF/B-VII evaluated discrete 
levels, and is shown in Fig. [TJtogether with the data from 
the present measurement. The maximum cross section 
for this level is roughly 60% higher in the ENDF/B-VII 
evaluation than in the present experimental data. The 
7 production cross section for this transition is shown in 
Fig. [15] The typical energy resolution for CdZnTe is 
about 2% at the Cd endpoint energy. In such a large en- 
ergy range, we can expect to see an additional 11, 12, and 
5 levels producing potential background 7 rays in the full 
energy peak (FEP), SEP, and DEP regions, respectively. 

e. 130 Te Endpoint Energy The endpoint of the 130 Te 
decay coincides with a ground-state transition from the 
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FIG. 15. The 7 production cross section for 2808-keV 7- 
ray in 63 Cu. This 7 ray transition is in the vicinity of the 
116 Cd endpoint energy. The corresponding level excitation 
cross section for this transition is shown in Fig. [7] and is 
included in the ENDF/B-VII evaluation for 63 Cu. 



12 



2533-keV level in 63 Cu; a level with significant strength 
(see Fig. [6]). The 7 production cross section for this 
transitions is shown in Fig. [16] This transition had the 
strongest 7 production cross section of any other tran- 
sition in this energy region. It is the strongest observed 
transitions above 1 MeV. Current double beta decay mea- 
surements involving Te typically employ TeC>2 crystals 
used as bolometers. There don't appear to be any ad- 
ditional troublesome 7-ray transitions from excited lev- 
els for this isotope, even when the energy window is ex- 
panded out to the anticipated 5-keV energy window. The 
2533-keV level has eight decay paths, with the dominate 
path being to the 1412-keV level in 63 Cu. For a dis- 
cussion of potential backgrounds from the 130 Te crystal 
please see Ref. [5B] . 




5 10 15 

Neutron Energy [MeV] 



FIG. 16. The 7 production cross section for 2535-keV 7- 
ray in 63 Cu. This 7 ray transition is in the vicinity of the 
130 Te endpoint energy. The corresponding level excitation 
cross section for this transition is shown in Fig. [6] and is 
included in the ENDF/B-VII evaluation for 63 Cu. 



/. 136 Xe Endpoint Energy Recently the EXO collab- 
oration reported the first ever observation of 2v/3f3 in 
136 Xe[57 . This measurement found that the endpoint of 
the decay to be at 2457.83(37) keV. There are no known 
7-ray transitions from excited levels in the immediate 
vicinity of the 136 Xe endpoint energy. The 3429-keV level 
has a 2468- keV 7 ray which is 11 keV away from the end- 
point of 136 Xe. This level primarily decays to the ground 
state, and 23% of the time decays to the second-excited 
level [42] . The level cross section for the 3429-keV level in 
63 Cu is shown if Fig. 17 
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FIG. 17. The 7 production cross section for 3429-keV 7-ray 
in 63 Cu. This 7 ray transition is in the vicinity of the 136 Xe 
endpoint energy, and is a ground state transition from the 
3429-keV energy level in 63 Cu. The corresponding level for 
this 7-ray is not included in the ENDF/B-VII evaluation for 
~Cu. 



interfering 7-ray lines within 10 keV of the energy of the 
single-escape peak. There are two 7-ray lines that are 
in 10-keV energy window corresponding to the energy of 
double escape peak. The 7 production cross section for 
the 3476-keV level in 63 Cu is shown in Fig. Ft8] This level 



Expanding the search window 
to include a 1.6% energy resolution reveals an additional 
3 levels with interfering 7-ray transitions. There are no 



only decays to the ground state, and is not included in 
ENDF/B-VII evaluation. 



g. 150 Nd Endpoint Energy There are a variety of ex- 
periments that are studying neutrinoless double beta de- 
cay in 150 Nd. The current experimental techniques to 
study this isotope rely on kTons of liquid scintillator, and 
estimate a 3% energy resolution 58 . Our analysis did 
not reveal any known levels with 7-ray transitions in the 
immediate vicinity of the endpoint energy of this decay. 
Table [I] gives the sensitivity of a 4-keV region around the 
3367-keV, or the endpoint for this decay. If there was a 
7 ray in this energy region, its maximum 7 production 
cross section would be below 0.30 mb. Expanding the 
search region to 3% energy resolution reveals an addi- 
tional 23 levels in nat Cu that produce full-energy 7 rays 
in the Nd region of interest. Again, as in the case of cal- 
cium, all of these levels are above the cutoff energy for 
the ENDF/B-VII evaluated discrete levels in either Cu 
nucleus. 
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FIG. 18. The 7 production cross section for 3476-keV 7-ray 
in 63 Cu. This 7 ray transition is in the vicinity of the 136 Xe 
double-escape peak energy region, and is a ground state tran- 
sition from the 3476-keV energy level in 63 Cu. The corre- 
sponding level for this 7-ray is not included in the ENDF/B- 
VII evaluation for 63 Cu 



IV. DISCUSSION AND CONCLUSIONS 

In conclusion, we have measured the 7 production 
cross sections for 111 transitions in nat Cu and eight addi- 
tional transitions in 62,64 Cu for neutron energies 1 MeV 
< E n < 100 MeV. We have compared our results with 
the ENDF/B-VII evaluations for Cu, and we found that 
the ENDF/B-VII evaluation agrees with our level cross 
sections for excited levels below 2 MeV for 63 Cu, and 
overestimates these level cross sections in 65 Cu. For en- 
ergy levels above 2 MeV we found significant discrepan- 
cies between the suggested level cross sections for both 
nuclei and our data. Our data were also compared to 
both the TALYS-1.2, and CoH3 statistical model calcula- 
tions. Both calculations displayed considerable variabil- 
ity in the agreement with the experimental data; achiev- 
ing excellent agreement for some levels, while differing 
by 300% for other excited levels. Given this degree of 
uncertainty in the models for most of the levels relevant 
for low-background, high-sensitivity studies need to be 
conducted to measure the level cross sections directly. 

In addition to examining discrete levels, we also com- 
pared the total neutron inelastic cross sections for both 



Cu nuclei with the ENDF/B-VII evaluations. We found 
excellent agreement for the 63 Cu total neutron inelas- 
tic cross section with the ENDF/B-VII evaluation and 
previous measurements. For the 65 Cu total neutron in- 
elastic cross section, our measurements were 30% lower 
than the ENDF/B-VII evaluation, which we attribute to 
unobserved transitions in 65 Cu. 

In addition to comparing our measurements with the 
ENDF/B-VII evaluation for Cu, we also looked at the im- 
plementation of these evaluations in the Monte Carlo pro- 
gram, Geant4. We found that the Geant4 was properly 
reproducing the overall shape and strength ENDF/B-VII 
evaluated cross section. On a more detailed level, how- 
ever we found that the decay properties of the nucleus 
were not being properly modeled to the degree neces- 
sary for estimating backgrounds in rare-event searches. 
There are relatively few discrete levels in ENDF/B-VII, 
and above 2 MeV the overwhelming strength of the con- 
tinuum begins to dominate the simulated inelastic inter- 
actions with the nucleus. In the current implementation, 
the ENDF /B-VII continuum file does not provide a decay 
path, nor do the de-excitation 7 rays correspond to tran- 
sitions in Cu. This approach needs to be reassessed for 
rare-decay searches, where the concern is 7-ray transi- 
tions that might interfere, or even replicate the signature 
of the experimental process. 

Finally, we examined the potential implications of our 
measurements on Oz//3/3 measurements. We found that 
of the nine frequently studied 0v(3/3 isotopes, only four 
had Q- values below the cutoff for ENDF/B-VII evalu- 
ated discrete levels in either Cu nucleus. We were able 
to identify 16 7-ray transitions in Cu as potential back- 
grounds for Ov/3/3 experiments. We measured level cross 
sections for nine of the levels that produce these 7-ray 
lines, and were able to put limits on the 7 production 
cross sections for regions where no specific 7-ray tran- 
sitions were identified. We identified several transitions 
that are problematic for experiments intending to use 
their 0i/fij3 isotope in a bolometer. 
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